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The transfer of turbulent energy through an inertial range from the driving scale to dissipative 
scales in a kinetic plasma followed by the conversion of this energy into heat is a fundamental plasma 
physics process. A theoretical foundation for the study of this process is constructed, but the details 
of the kinetic cascade are not well understood. Several important properties are identified: (a) the 
conservation of a generalized energy by the cascade; (b) the need for collisions to increase entropy 
and realize irreversible plasma heating; and (c) the key role played by the entropy cascade — a dual 
cascade of energy to small scales in both physical and velocity space — to convert ultimately the 
turbulent energy into heat. A strategy for nonlinear numerical simulations of kinetic turbulence is 
outlined. Initial numerical results are consistent with the operation of the entropy cascade. Inertial 
range turbulence arises in a broad range of space and astrophysical plasmas and may play an 
important role in the thermalization of fusion energy in burning plasmas. 
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I. INTRODUCTION 

In a warm, magnetized kinetic plasma, the dissipation 
of low-frequency turbulent fluctuations by kinetic mech- 
anisms occurs most strongly at length scales of order the 
plasma particle Larmor radii. If turbulence is driven in 
the plasma at a much larger scale, an inertial range devel- 
ops in which the injected energy is cascaded by nonlinear 
interactions to the small scales at which the turbulence 
can be effectively dissipated. This transfer of turbulent 
energy through an inertial range from the driving scale 
to dissipative scales followed by the conversion of this 
energy into heat is a fundamental kinetic plasma physics 
process. In this paper, I lay out a theoretical foundation 
on which to build our knowledge of inertial range turbu- 
lence in kinetic plasmas and outline a practical strategy 
for the nonlinear numerical modeling of this process. 

The study of turbulence in kinetic plasmas is not a new 
endeavor; in magnetic fusion community, understanding 
the effect of turbulence on the transport properties of the 
plasma is one of the primary goals of the numerical mod- 
eling effort. But the turbulence most often investigated 
in laboratory plasmas is driven by microinstabilities aris- 
ing from gradients in the background plasma — for exam- 
ple, ion temperature gradient (ITG) and electron teinper- 
ature gradient (ETG) instabilities [H, i, i, i, i, i, 0, i, H ■ 
These instabilities inject energy into the plasma at the 
scale of the particle Larmor radius; the resulting tur- 
bulence does not develop an inertial range because the 
driving and kinetic dissipation scales are of roughly the 
same order. Inertial range turbulence in kinetic plas- 
mas, therefore, represents a fundamentally different pro- 
cess than the kinetic turbulence traditionally studied by 
the fusion community. 

The turbulence in most astrophysical contexts, on the 
other hand, is typically driven at scales much larger than 
the Larmor radius and gives rise to a turbulent cas- 
cade through an inertial range to the dissipative scales at 



which the energy of turbulent ffuctuations is ultimately 
converted into particle thermal energy. Before construct- 
ing the detailed theoretical foundation for the study of 
this process, presented in Sec |TT1 I will provide here a 
simplified blueprint of the physical mechanisms guiding 
the flow of energy. It is important to emphasize that 
many of the mechanisms described below are not well 
understood. Here I aim only to paint the overall pic- 
ture in broad strokes, deferring discussion of the many 
uncertainties to SecHTl 

At a scale larger than the particle mean free path — 
a collisional scale for which magnetohydrodynamics 
(MHD) provides an adequate description — the plasma 
is stirred by some external mechanism, driving an as- 
sortment of MHD Alfven, fast, slow, and entropy mode 
fluctuations in the plasma. At this driving scale, these 
modes are undamped; thus, a turbulent cascade develops 
nonlinearly to transfer the fluctuation energy to smaller 
scales. The compressive modes become damped as the 
cascade reaches scales of order or smaller than the col- 
lisional mean free path, but the Alfvenic cascade con- 
tinues undamped down to the scale of the ion Larmor 
radius. At this kinetic scale, the electromagnetic fluc- 
tuations may be damped coUisionlessly by the Landau 
resonance with the ions. In the absence of collisions, this 
process conserves a generalized energy — the free energy 
removed from the electromagnetic fluctuations generates 
nonthermal structure in velocity space of the ion distri- 
bution function. The remaining electromagnetic fluctua- 
tion energy continues to cascade below the scale of the ion 
Larmor radius as a kinetic Alfven wave cascade. Upon 
reaching the scale of the electron Larmor radius, the elec- 
tromagnetic fluctuations of the kinetic Alfven wave cas- 
cade are completely damped via the Landau resonance 
with the electrons; again, a generalized energy is con- 
served in this process, leading to the creation of nonther- 
mal structure in velocity space of the electrons. But the 
damping of the electromagnetic fluctuations and conse- 
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quent generation of structure in velocity space does not 
correspond to heating; irreversible heating requires an in- 
crease in entropy that can only be achieved by collisions. 
The thermalization of the turbulent energy by collisions 
is ultimately achieved thanks to a cascade to small scales 
in velocity space of the particle distribution functions — 
an entropy cascade. The entropy cascade drives the dis- 
tribution function structure in velocity space to scales 
small enough that even weak collisions are sufficient to 
smooth out that structure towards the Maxwellian, caus- 
ing entropy to increase — this is the final step in the con- 
version of the energy of the turbulent fluctuations to ther- 
mal energy of the plasma particles. 

The entire process described above is the kinetic gener- 
alization of the familiar cascade of energy in a fluid turbu- 
lent system; this fundamental kinetic plasma physics pro- 
cess, encompassing the dual cascade in both physical and 
velocity space, is referred to as the kinetic cascade of the 
generalized energy [l^] ■ Neither the detailed interactions 
of the kinetic cascade nor its implications for the flow of 
energy in turbulent systems are well understood; a theo- 
retical foundation for its study is layed out in SecHIl Non- 
linear numerical simulations of the kinetic cascade occur- 
ring in inertial range turbulence are expected to play a 
leading role in shedding light on the complex physical 
mechanisms involved — a strategy for the numerical ap- 
proach to this problem is presented in Sec IIIII and initial 
numerical results are discussed in Sec IIVI Astrophysical 
and laboratory environments in which the kinetic cascade 
of inertial range turbulence may play an important role 
are identified in SecfVl 

II. THEORETICAL FOUNDATION 

Although inertial range turbulence plays an important 
role in the evolution of many kinetic plasmas, it has not 
yet attracted great attention within the scientific commu- 
nity. The physics underlying this intersection of kinetic 
plasmas and turbulence is fascinating yet remains poorly 
understood. This paper aims to construct a foundation 
upon which a thorough knowledge of the varied mech- 
anisms involved may be built. I present here a view of 
the overall context in which the kinetic cascade operates 
and aim to emphasize the uncertainty underlying some 
aspects of this view. 

I begin with the assumption that the turbulence is 
driven at a scale L much larger than the ion coUisional 
mean free path, L 3> ^nifpr "^^^ dynamics at the driv- 
ing scale is therefore collisionally dominated, so the single 
fiuid magnetohydrodynamic (MHD) theory provides an 
adequate description of the plasma dynamics at all par- 
allel scales /|| satisfying ^|| ^ -^mfpi (^here _L and || in 
this paper are relative to the direction of the local mean 
magnetic field) [s^l . The assumption L ^ ^uitpi ^^^^ '^'-'^ 
be true for all environments in which inertial range, ki- 
netic turbulence arises; the implications of relaxing this 
assumption are discussed in Sec III II Based on the fact 



that the damping of all linear wave modes at the MHD 
scales Z|| ^ '^mfpi negligible, the damping of the tur- 
bulent fluctuations is likewise expected to be negligible; 
an inertial range will then develop to cascade the tur- 
bulent energy via nonlinear interactions to ever smaller 
scales. This cascade will continue down to the transi- 
tion from fluid to kinetic behavior at a scale /y ~ '^mfpr 
At this point, a fluid treatment of the turbulence is no 
longer adequate — instead a kinetic theory must be used 
to describe the dynamics. The nature of the kinetic con- 
tinuation of the inertial range to scales Z|| < Aj^^p^ de- 
pends critically on the character of the MHD turbulent 
fluctuations reaching this transition. 

Here I review five key concepts underpinning the the- 
ory of MHD turbulence before applying them to con- 
struct a model for the fluid portion of the turbulent in- 
ertial range. 

1. The Kolmogorov Hypothesis: the nonlinear energy 
transfer is constant through all scales within the 
inertial range, defined as the extent of scales influ- 
enced neither by the energy injection mechanism 
nor by dissipation; at a given scale, the energy 
transfer is spectrally local, with a rate determined 
only by the turbulent conditions at that scale 

2. The Kraichnan Hypothesis: even in a magnetized 
plasma with no mean field, the magnetic field 
of large-scale fluctuations behaves effectively as a 
mean fleld for fluctuations on smaller scales p^ . 

3. Colliding Wavepackets: as is apparent when the 
equations of incompressible MHD are cast into 
Elsasser form, only Alfven wave packets traveling in 
opposite directions along the mean magnetic field 
interact nonlinearly 

4. Critical Balance: as strong MHD turbulence cas- 
cades to small scale, it maintains a state of balance 
between the (parallel) linear propagation and (per- 
pendicular) nonlinear interaction timescales [3, 

5. Decoupling of Fast Wave Dynamics: theory and 
simulations of compressible MHD turbulence sug- 
gest that the isotropic cascade of fast MHD wave 
turbulence decouples from the dynamics of the 
Alfven, slow, and entropy modes [16. ilTij. 

A. Weak MHD Turbulence 

Adopting the Kolmogorov and Kraichnan hypotheses, 
the interaction of oppositely directed, small-amplitude 
Alfven wave packets along the mean (or large-scale) mag- 
netic field in incompressible MHD is considered. If it 
is assumed (a) that the nonlinear interaction between 
two colliding wavepackets is weak — and therefore that a 
single wave packet must undergo many collisions before 
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it has been distorted enough to have effectively trans- 
ferred its energy to smaller scale — and (b) that the energy 
is transferred to smaller scales isotropically, an inertial 
range one-dimensional energy spectrum that scales with 
wavenumber as fc^^/^ is predicted. A spectrum with this 
scaling is often referred to as the Iroshnikov-Kraichnan 
spectrum [l^, • 

However, evidence from laboratory plasmas [16 
l20l |. the solar wind [2l[ and numerical simulations [22 
suggests that the MHD turbulent cascade is not isotropic 
but preferentially transfers energy to small scales per- 
pendicular to the mean magnetic field. Taking the 
observed anisotropy into account, anisotropic theories 
of MHD turbulence have been proposed and refined 
[il HI, m m, iS [13, Hi] , and have led to the emer- 
gence of a mature theory for weak incompressible MHD 
turbulence [1^ [13, [3l| . The predicted one-dimensional 
energy spectrum for weak incompressible MHD turbu- 
lence scales as fcj^; the cascade is anisotropic, with 
all energy cascading to higher perpendicular wavenum- 
bers and none to higher parallel wavenumbers. Un- 
like the isotropic Iroshnikov-Kraichnan cascade, in which 
the nonlinear interactions weaken as the turbulence pro- 
gresses to smaller scale, this anisotropic theory predicts 
that, as the turbulence cascades to higher perpendicu- 
lar wavenumber, the nonlinear interactions strengthen, 
leading eventually to the inevitable violation of the as- 
sumption of weak nonlinear interactions [1^ [H, [sO, [IH . 
Thus arises the most important implication of the weak 
turbulence theory: given a cascade through a sufhciently 
broad range of scales, weak incompressible MHD turbu- 
lence will inevitably transition to a state of strong tur- 
bulence. 



B. Strong MHD Turbulence 

Below the perpendicular scale where the perturbative 
treatment fails, the resulting state of strong incompress- 
ible MHD turbulence is treated phenomenologically by 
assuming that the turbulence maintains a state of criti- 
cal balance, lu ~ tOnh as it cascades to smaller scales. For 
a linear Alfven wave frequency, lo = k^^VA, and a nonlin- 
ear frequency determined by the perpendicular dynamics, 
w„i ~ k±v±, theory predicts a one-dimensional kinetic 
energy spectrum Ek{kj_) cx k^^^^ and a scale-dependent 

anisotropy fcy oc fc^^'^. Although numerical simulations of 
strong MHD turbulence with a dynamically strong mean 
field appear to to support the scale-dependent anisotropy 
k\\ cx A:^^'^ [33,[3l|, these same simulations routinely pro- 
duce one-dimensional energy spectra that appear to scale 
as fcj^^^^ rather than k^^^^ [s^, [s^. A modification of 
the theory to take into account a proposed polarization 
alignment f35| may shed light on this discrepancy, but 
the matter of strong MHD turbulence remains controver- 
sial. It is also worthwhile noting here that the predicted 
transition from weak to strong MHD turbulence has only 



recently been numerically reproduced \?>u\ . 

In spite of the uncertainties in the theory for the strong 
incompressible MHD turbulence, the combination of key 
concepts (l)-(4) in Sec HI] leads to the following general 
prediction: in any magnetized, incompressible plasma, 
whether or not a strong mean field exists at the large 
scale, at sufficiently small scales a state of strong MHD 
turbulence will arise. Therefore, given the assumption of 
sufficiently large driving scale L ^ -^mfpi' transi- 
tion from fluid to kinetic behavior l\\ ~ '^nifpi' turbu- 
lent fluctuations will be dominantly perpendicular with 
/c|| ^ fcj^ and the fluctuation amplitudes will be small 
compared to the mean field, 5B±^ ^ Bq. The kinetic 
description of the continued turbulent cascade to smaller 
scales is strongly influenced by these predicted charac- 
teristics of the turbulence. 



C. Compressible MHD Turbulence 

The discussion thus far has been limited to incom- 
pressible MHD. In the case of compressible MHD, an 
arbitrary stirring mechanism may inject energy into the 
compressive modes of the plasma — the fast, slow, and 
entropy modes. A weak turbulence treatment of com- 
pressible MHD for fast and Alfven waves suggests that 
only a small amount of energy is transferred from fast 
waves to Alfven waves at large A;|| [sj. Numerical simula- 
tions of strong turbulence in compressible MHD demon- 
strate an isotropic cascade of fast waves that scales as 
/c^^/^ and suggest that this cascade of fast wave en- 
ergy is decoupled from the Alfven and slow wave cas- 
cades • In the region of wavevector space where the 
the energy of the Alfven wave cascade is conc entrated , 
k\\ <C k±^, the fast wave frequencies uj ^ k^y^c^ + ' 
where = jp/iriiUii) is the sound speed, are much 
higher than the Alfven wave frequencies uj = k\\VA, so 
the decoupling of these cascades is not surprising. The 
cascade of the slow waves and entropy modes is slaved 
to the Alfven wave cascade [l^ [s^ and is seen to mimic 
the anisotropic cascade of the Alfven waves [13, [1^ . Al- 
though, from these arguments, the compressive modes of 
the plasma turbulence are not expected to significantly 
alter the behavior of the Alfven wave cascade as deter- 
mined by incompressible MHD turbulence theory, this 
topic remains an open area of research. 



D. Transition to Kinetic Turbulence at l\\ 



^mfpi 



As the inertial range turbulent cascade enters the 
regime of weak collisionality at /j| ^ -^mfpi' nature 
of the fluctuations at this transition motivates a partic- 
ular choice for the kinetic description of the dynamics of 
the cascade to yet smaller scales. The theory of MHD 
turbulence described above predicts the following turbu- 
lent state at the transition: isotropic fast wave fluctu- 
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ations and anisotropic Alfven, slow, and entropy mode 
fluctuations with fe|| <gC k±. Because both the energy of 
the fluctuations decreases with increasing wavenumber 
(oc for fast waves and oc kj_^^^ for Alfven, slow, 

and entropy modes) and the driving scale is assumed 
large L ^ -^rnfpi' amplitude of the magnetic field 
fluctuations at this scale will be small compared to the 
local mean field, 5B± <C -Bo- Although strong nonlin- 
ear interactions dominate the turbulence — leading to the 
nonlinear cascade of the energy of a wave on a timescale 
of order the wave period — the small amplitude of the fluc- 
tuations suggests that it is not unreasonable to estimate 
the coUisionless damping of the turbulent fluctuations us- 
ing the damping rates from linear theory. For plasmas 
with (3i > 1, linear theory shows that the kinetic wave 
modes corresponding to the fast and slow MHD wave 
modes are significantly damped coUisionlessly at scales 
^11 ^ '^mfpi IH) Soli for lower (3i plasmas, however, 
the damping is less vigorous. Fast waves are also subject 
to strong dissipation as they steepen into shocks. Al- 
though the argument that these wave modes are damped 
rapidly in the regime of weak coUisionality is plausible, 
further exploration of this matter is required. For the 
remainder of this paper, I will neglect the effect of any 
energy in the fast wave mode (and its kinetic continua- 
tion) on the dynamics of the Alfven, entropy, and slow 
modes. 

In this limit k\i ^ k± and SB± ^ Bq, the kinetic tur- 
bulent dynamics of the Alfven, slow, and entropy modes 
is rigorously described by a low-frequenc y e xpansion of 
kinetic theory called gyrokinetics [13, Hi], [H, 123, HI, HI, 
|46|, |47|, |48|, |49|, |50|, IHf . Exploiting the timescale separa- 
tion between frequency of the turbulent fluctuations and 
the ion cyclotron frequency, w <C f^i, gyrokinetics aver- 
ages over the fast cyclotron motion of charged particles 
in the mean magnetic field. The gyrokinetic approxima- 
tion orders out the fast MHD wave and the cyclotron 
resonance, but retains finite-Larmor-radius effects, col- 
lisionless damping via the Landau resonance, and colli- 
sions. A powerful result from gyrokinetic theory is the 
conservation of a generalized energy W according to the 
equation 
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where the index s denotes plasma species, 5fs is the per- 
turbation from the Maxwellian equilibrium distribution 
function Fqs with temperature Tqs, hg is the gyrocenter 
distribution function defined by 5fs — —qs4>Fos/Tos + hs 
with (j) as the scalar potential, 3a is an external antenna 
current drivin g th e system, and (. . .)c denotes the colli- 
sion operator [l^, |50, [sj . This relation has important 
implications for the flow of energy in kinetic turbulence, 
as discussed in detail below. 



E. Kinetic Turbulence at k±pi <C 1 



The dynamics of kinetic turbulence for all scales Z|| 



< 



'^mfpi described by gyrokinetics as long as the fre- 
quency of fluctuations remains below the ion cyclotron 
frequency [5^ ; the failure of the gyrokinetic theory when 
w ^ i7i is discussed in Sec IH II Here I consider the prop- 
erties of this turbulence in the kinetic regime at perpen- 
dicular scales larger than the ion Larmor radius, l± 
or k±pi 1. In this regime, Schekochihin et al. [lO| 
have demonstrated that the turbulence has the follow- 
ing characteristics: (a) the dynamics of the Alfven wave 
cascade decouple from the compressive slow wave and en- 
tropy mode fluctuations and ar e rig orously described by 
the equations of reduced MHD [2a[2i.[54|: (b) the slow 
and entropy mode dynamics is governed by an ion ki- 
netic equation and these modes are passively mixed non- 
linearly by the Alfven waves; and (c) no energy is trans- 
ferred between the Alfvenic fluctuations and the com- 
pressive fluctuations. Since the Alfven wave dynamics 
are given by reduced MHD, the Alfven waves undergo 
the same anisotropic cascade described by the theory for 
strong incompressible MHD turbulence and the cascade 
continues undamped down to the perpendicular scale of 
the ion Larmor radius k^^pi ^ 1. Whether the slow and 
entropy fluctuations are damped in this regime or merely 
passive mixed by the Alfven waves with little damping is 
uncertain at present [IQ] . 



F. Kinetic Turbulence at k\ 



pi 



At the perpendicular scale of the ion Larmor radius 
k±Pi ^ 1, two effects come into play. First, the electro- 
magnetic fluctuations may be damped via the Landau 
resonance with the ions. In the absence of collisions or an 
external driving mechanism, equation ^ demonstrates 
that the generalized energy W must be conserved dur- 
ing this process. Through the wave-particle interaction, 
the electromagnetic fluctuations cause parallel accelera- 
tion of the ions, transferring energy from the fields to 
the ions — the energy lost by the electromagnetic fields is 
converted into nonthermal structure in velocity space of 
the ion distribution function. 

The second effect is the decoupling of the ion motion 
from the electromagnetic fluctuations. An ion samples 
the electromagnetic field over the scale swept out by its 
fast Larmor motion. For field fluctuations on length 
scales smaller than the ion Larmor radius k^pi > 1, 
this leads to an averaging over the spatially oscillating 
field; the net field experienced by the particle decreases 
through the ring-averaging over the Larmor motion ^S^l . 

This gradual decoupling of the ions from the fields 
leads to an inherently nonlinear effect — the entropy cas- 
cade. Only recently identified [l0[, the entropy cascade 
plays a critical role in the inevitable thermalization of the 
turbulent energy. Due to the fact that the radius of the 
Larmor motion for an individual particle depends on the 
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perpendicular kinetic velocity v± of that particle, parti- 
cles with high v± decouple from the field more rapidly 
than those with low v± . The drift velocity of each par- 
ticle perpendicular to the magnetic field is the E x B 
velocity; particles with low perpendicular kinetic veloc- 
ities v± feel a stronger ring-averaged electric field (E) 
than particles with higher and so drift faster. This 
differential (E) x B drift leads to a coupling of structure 
in physical space with that in velocity space, resulting in 
nonlinear phase mixing to smaller scales in both physi- 
cal and velocity space. This dual cascade is referred to 
as the entropy cascade and occurs for any species s at 
scales k±pi < 1. As will be explained in Sec III HI this 
inherently nonlinear phenomenon plays a critical role in 
the thermalization of turbulent energy. 



where ^ = v\\/v is the pitch angle and is a colli- 
sional frequency. For weakly coUisional plasmas, there- 
fore, the nonthermal structure in velocity space of hg 
(created through the parallel acceleration of particles by 
electromagnetic waves via the Landau resonance) must 
be driven to sufficiently small scales in velocity space such 
that vd^{hs/FQs)/d^^ ~ lo; only then can collisions act 
rapidly enough to diffuse the fine scale structure in veloc- 
ity space and increase the entropy. The entropy cascade 
is therefore the mechanism by which velocity-space struc- 
ture is driven to very small scales at which point the free 
energy contained in the nonthermal fluctuations can be 
thermalized. This ultimate thermalization of the turbu- 
lent energy marks the endpoint of the kinetic cascade. 



G. Kinetic Turbulence at k±pi S> 1 

Any energy remaining in the electromagnetic fluctua- 
tions continues to cascade to scales smaller than the ion 
Larmor radius fc_LPi ^ 1 as a kinetic Alfven wave cas- 
cade. For many parameter regimes, the damping of the 
turbulent fluctuations via the Landau resonance with the 
electrons is often substantial for the entire range of scales 
k^Pi ^ 1 .53] • Regardless of the parameter regime, how- 
ever, by the time the cascade reaches the scale of the 
electron Larmor radius, /c^pe ~ 1, the damping always 
becomes strong enough to terminate the cascade via the 
Landau resonance with the electrons. Again, in the ab- 
sence of collisions, the generalized energy W is conserved 
and the energy lost by the fields builds nonthermal struc- 
ture in velocity space of the electron distribution func- 
tion. 



H. Thermalization of the Turbulent Energy 

As demonstrated by Howes et al. [50] , damping of the 
electromagnetic fluctuations by the Landau resonance 
does not represent heating of the plasma; collisions are 
required to increase the entropy of the system and real- 
ize irreversible plasma heating. The change in entropy 
Ss can be expressed as 

where V is the plasma volume. It can be shown that the 
integral is non-positive, so the entropy always increases. 
Comparing to equation it is clear that, in the absence 
of driving, the generalized energy may only decrease due 
to collisions; this collision term represents an increase in 
entropy and thus represents irreversible heating of the 
plasma. 

The action of the collision operator is to smooth out 
structure in velocity space of the distribution function; 
the collision operator is proportional to vd'^{hs/FQs)/d^'^ 



I. Other Considerations 

Gyrokinetic theory is only valid for dynamics with 
characteristic fluctuation frequencies small compared to 
the ion cyclotron frequency, uj <^ Vli. Since the fluctua- 
tion frequency typically increases as the turbulence cas- 
cades to smaller scales, it is possible to violate this condi- 
tion at a scale deep into the kinetic cascade at wavenum- 
bers k^pi 3> 1; however, coUisionless damping of the 
turbulent fluctuations can slow this frequency increase 
along the cascade [53]. Neglecting this damping provides 
a conservative estimate of the perpendicular wavenumber 
threshold at which the cyclotron frequency is reached; in 
this case, the frequency can be estimated throughout the 
kinetic regime by 

/ / \ 2/3 

1^ VA + 2/(i + Te/r,); ^l) 

[531 ] . The small scale limit of applicability of gyrokinetic 
theory must be determined on a case by case basis — see 
Howes et al. [s^ for a detailed discussion of the limits of 
validity of gyrokinetics in the solar wind. Astrophysical 
plasmas in general support an inertial range of many or- 
ders or magnitude, so this condition is usually satisfied to 
wavenumbers beyond the ultimate termination of the ki- 
netic cascade at fcj^pe ~ 1- Taking the specific example of 
a plasma with = 1 and Ti/Tc. = 1, this low- frequency 
requirement is met when ki_pi <C (Lj pi)^^'^\ for a typi- 
cal astrophysical inertial range of Lj pi ~ 10^, this gives 
k^Pi <^ 100, so the ion cyclotron frequency is not reached 
until beyond the scale of k±pe ~ 1 (equal to k±pi ~ 40). 

Next I consider how this picture of the kinetic cascade 
changes if the cascade is driven at a coUisionless scale 
L < Ajjjfp^. Since this work is concerned with inertial 
range turbulence, it is implicitly assumed that the driv- 
ing scale is much larger than the kinetic scales, L ^ pi. 
Based on the fact that the Alfven wave cascade, even 
in coUisionless regimes is governed by fluid equations for 
all scales k±pi ^ 1, it seems reasonable to expect that 
the predictions of strong incompressible MHD turbulence 
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would continue to hold: an anisotropic cascade will de- 
velop leading to an energy spectrum that scales as k^^^ 

and a scale-dependent anisotropy k\\ cx /c^^'^. If the tur- 
bulence is driven isotropically at scale L, one expects that 
gyrokinetic approximation will be satisfied for all scales 
such that {k^LY^^ ^ 1 — in other words, the specific 
predictions for the kinetic cascade outlined here, based 
on the gyrokinetic theory, will apply beginning at a scale 
several orders of magnitude below the driving scale. 

III. NUMERICAL APPROACH 

The efficient numerical study of inertial range turbu- 
lence in kinetic plasmas is best conducted using a mixture 
of computational approaches. At the large, coUisional 
scales, a fluid theory is suitable; a hybrid scheme of fluid 
electrons and kinetic ions is most efficient at intermedi- 
ate scales; at the small scales where the turbulent fluc- 
tuations are damped and their energy is converted into 
heat, a kinetic scheme is required for both species. A 
hierarchy of such schemes is rigorously derived from ki- 
netic theory by Schckochihin et al. [l^. In this section, 
I describe a numerical strategy for the study of the most 
fascinating aspect of inertial range turbulence in kinetic 
plasmas — the transition regime from the (reduced) MHD 
Alfven wave cascade to the kinetic Alfven wave cascade 
at the scale of the ion Larmor radius. 

The need for a kinetic treatment of the transition 
regime limits the spatial dynamic range possible in a nu- 
merical simulation. The strategy chosen is to model both 
ions and electrons gyrokinetically in the neighborhood of 
k^Pi ~ 1. The simulation results reported here were 
obtained using AstroGK, a new gyrokinetic simulation 
code designed for the numerical modeling of astrophys- 
ical turbulence. The code is based on GS2, a mature, 
widely used gyrokinetic code for the design and inter- 
pretation of laboratory experiments in the magnetic fu- 
sion program 0, [s^. AstroGK is an Eulerian flux tube 
code with periodic boundary conditions that evolves the 
five-dimensional distribution function for each plasma 
species. The spatial components x and y are handled 
spectrally in Fourier space, and the z-direction is handled 
with a compact finite differencing scheme. Integration 
over the velocity-space grids is accomplished using spec- 
tral integration by quadrature in energy and pitch an- 
gle, while the differentiation needed by the momentum- 
conserving collision operator in pitch angle is accom- 
plished using finite differencing. The positions of the 
energy and pitch angle grid points are determined us- 
ing Legendre polynomials — an example with 128 points 
in pitch angle and 32 points in energy is shown in Fig- 
ure [T] The dynamics are fully electromagnetic, evolv- 
ing the scalar potential (p, the parallel component of the 
vector potential Ay , and the parallel component of per- 
turbed magnetic field 5B\\ to describe the fiuctuating 
electromagnetic fields. The linear terms are advanced 
implicitly to avoid the need to satisfy a Courant con- 
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FIG. 1: Plot of grid used in velocity space with 128 pitch 
angles and 32 energies. Grid point locations are chosen spec- 
trally in a Legendre polynomial basis. 



dition for the fast electron dynamics. Nonlinear terms 
are evaluated by fast Fourier transform in real space and 
the term is advanced using a 3rd-order Adams-Bashforth 
scheme. A detailed description of the algorithms used 
within AstroGK is given in Howes et al. [56|. 

Because the range of scales comprising the transition 
regime falls within the larger kinetic cascade, two devel- 
opments are essential for the successful modeling of this 
system: (a) a scheme for modeling the inertial range en- 
ergy transfer from turbulent fiuctuations at scales larger 
than the simulation domain; and (b) a mechanism for 
dissipating the turbulence artificially to avoid the un- 
physical buildup of energy at the smallest resolved scales. 
Both of these developments will be described in detail 
and validated in a future work, but here we provide a 
brief description. 

Critical balance posits that the timescale for the non- 
linear transfer of energy through the inertial range is of 
order the linear wave frequency at that scale in the cas- 
cade (see Sec HI]). To simulate this process, the plasma 
is driven by an external antenna at several large wave- 
lengths comparable to the scale of simulation domain — it 
is important to excite separately traveling waves moving 
in both directions along the mean field in order to realize 
the colliding wave packets necessary for turbulence. The 
amplitude and phase of each mode determined accord- 
ing to a Langevin equation with a decorrelation time of 
order the wave period to emulate the conditions within 
the inertial range of strong MHD turbulence. The an- 
tenna drives only a perpendicular fiuctuating magnetic 
field to excite primarily Alfven waves (although at the 
amplitudes needed to drive strong turbulence, fiuctua- 
tions in the direction of the local magnetic field lead to 
some power coupled into compressive fluctuations). 

Although wave-particle interaction through the Lan- 
dau resonance acts to damp the turbulent cascade, it is 
often insufficient to terminate the cascade within the re- 
solved range of scales. To prevent a bottleneck of en- 
ergy at the smallest scales, we act on the turbulence 
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with a hypercoUisional operator — this operator acts on 
the distribution function with a non-momentum conserv- 
ing coUision operator using a scale-dependent coefficient 
i'ffs(fcj_/fc±max)f , where ps S [4,8]. The heating caused 
by hypercolhsionahty is positive definite; with a carefully 
chosen coefficient VHe^ the electron hypercoUisionality is 
sufficient to terminate the cascade on resolved scales. 

Choosing the ideal value of uue for a given simu- 
lation is often difficult. To avoid resorting to trail 
and error, an adaptive scheme is used to choose and 
modify the electron hypercoUisional coefficient based on 
nonlinear estimations of the energy transfer frequency 
and of the effective coUisional (including the hypercol- 
lisional term) damping rate at a given scale k±HPi- 
The energy transfer frequency at k±HPi is estimated by 
Wni ^ k±H^ \ SB±{kx, fcj,)p/87r where the sum is over all 
Fourier modes (kx,ky) within a band of width Ak±HPi 
centered at kj^nPi- The nonlinear damping rate is given 
by the total coUisional and hypercoUisional heating [50] 
within the same perpendicular wavenumber band divided 
by the sum of the total fluctuation energy in the wave 

band, 7„, ~ J2(.Pc^ + Pee + PHcr + PHce)/T.E- The 

coefficient of electron hypercoUisionality v^e is then ad- 
justed within preset bounds in order to achieve a value 
oi "fni/tOni — l/(27r) within some tolerance. 

Motivated by the desire to mitigate the large compu- 
tational cost of gyrokinetic simulations of kinetic turbu- 
lence, a scheme for reducing the time to reach a statisti- 
cally steady state of turbulence is used. Initially, a small 
simulation is run for several eddy turn-around times at 
the driving scale. The number of cpu-hours required to 
achieve a steady-state for the small problem is moder- 
ate. The restart files from this initial run are processed 
to add more spatial Fourier modes in the perpendicular 
direction and finite-difference grid points in the parallel 
direction — this effectively increases the spatial resolution 
of the simulation. The initial values of the added Fourier 
modes are zero and the values at the additional parallel 
grid points are determined by spline interpolation. When 
the simulation is restarted, to reach a new steady state, 
it only needs to run for an eddy turnaround time of the 
smallest non-zero mode. This Fourier mode expansion 
technique is used recursively until the desired resolution 
is reached. 



IV. NUMERICAL RESULTS 

Here we present results from a simulation of the tran- 
sition regime from the MHD Alfven wave cascade to the 
kinetic Alfven wave cascade for a plasma with — 1 and 
Ti/Tg — 1. The spatial dimensions in the plane perpen- 
dicular to the mean field are treated pseudospectrally 
on a 32 X 32 grid; the parallel direction is treated us- 
ing a compact finite-difference scheme on 64 grid points. 
Velocity space resolution uses 128 points in pitch angle 
and 32 points in energy, for a total of 4096 points over 
the half-plane uj^ > in velocity space; the positions 



of grid points are chosen pseudospectrally for maximum 
accuracy when performing integrations over velocity [56j . 
The distribution of velocity space grid points for this sim- 
ulation is shown in Figure [TJ A fully ionized proton and 
electron plasma is specified with a realistic mass ratio of 
mi/rrie = 1836 and both species are treated gyrokineti- 
cally. In summary, the dimensions of the simulation are 
(n^;, riy, riz, n^, n_E, «s) = (32, 32, 64, 128, 32, 2) for a total 
of 536,870,912 computational mesh points. The simula- 
tion used a total of 19,118 cpu-hours on Franklin, the 
Cray XT4 at NERSC. 

A momentum conserving pitch-angle scattering colli- 
sion operator is employed for like-species collisions with 
Vi = 0.001 and v,, — 0.001; interspecies collisions are 
neglected. The fully dealiased range of perpendicular 
wavenumbers covers fcj: Pi € [0.4,4.0] and fc^p; € [0.4,4.0]. 

The simulation is brought to the steady-state 
using the Fourier mode expansion technique de- 
scribed in Sec IIIIl the problem is run using a 
spatial grid of (16,16,32) for 1.3 outer scale peri- 
ods before expanding to (32,32,64). Zero initial 
conditions are specified. The Langevin antenna, as 
described in Sec IIIIl drives the six lowest wavenum- 
ber modes corresponding to (fc^j/fco, fcj,/fco, fcz/^zo) = 
(1, 0, 1), (0, 1, 1), (-1,0, 1), (1, 0,-1), (0, 1, -1), (-1,0,-1) 
with an amplitude of 5.0, a frequency oi lj = fc||W^ and 
a decorrelation frequency of Q.Sk\\VA- 

A fixed ion hypercoUisionality is used with vui = 0.04 
and Pi = 8. Electron hypercoUisionality is chosen adap- 
tively to terminate the cascade as described in Sec IIIIl 
In the initial (pre-expansion) run, the electron hyper- 
coUisionality is allowed to adjust within the bounds 
10 < i^He < 1200 to meet a dissipation criterion of 
0.16 < ^ni/uJni < 0.24 for aU modes 1.6 < k±pi < 2.4; 
in the expanded restart, the bounds and dissipation cri- 
terion are the same, but modes tested fall in the range 
3.6 < kj^pi < 4.4. The exponent for the electron hyper- 
coUisionality is = 8. 

The high velocity-space resolution simulation pre- 
sented here was set up to investigate the entropy cascade. 
Presented in Figure[l]are surface plots and projected con- 
tours over velocity space of the real part of the perturbed 
ion distribution function evolved by AstroGK, gi, defined 

by 

(4) 

Here the subscript k denotes a Fourier coefficient and 
{4>)-[isk is the gyro-average at the particle guiding cen- 
ter R^. The velocity space structure at the mid-plane 
in the parallel direction is shown for two Fourier modes 
{kxPi, kyPi). The upper plot shows the mode (0.4,0.4), 
corresponding to k±pi ~ 0.566; the lower plot (2.0,2.0), 
k±pi ~ 2.83. Wave-particle interactions via the Landau 
resonance drive structure in in velocity space; the non- 
linear phase mixing of the entropy cascade, on the other 
hand, creates structure in the v± direction. The mode 
at k±pi ~ 0.566 in the upper plot shows primarily struc- 
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FIG. 2: Surface plots and projected contours of the structure 
of the real part of the ion distribution function Qs in velocity 
space for the Fourier modes (kxPi,kypi) = (0.4,0.4) (upper 
panel) and (2.0, 2.0) (lower panel) at the mid-plane of the 
simulation domain in z. The lower plot, at a wavenumber 
k^pi > 1, shows v±_ structure in velocity-space, a finding 
consistent with the presence of an entropy cascade. 



ture in W|| with little variation in v^; the lower plot, at 
k^Pi ~ 2.83, shows more w_l structure in velocity space, a 
finding consistent with the action of the entropy cascade 
at k_ipi > 1. Velocity-space in the electron distribution 
function for these same Fourier modes (not shown) is 
much smoother, as expected for k^p^ <C 1. 

Although the results presented here are insufficient to 
identify clearly the ion entropy cascade, they represent a 
promising first step. The development of more sophisti- 
cated diagnostics of velocity space — for example, taking 
spectral transforms in velocity space to yield a velocity- 
space power spectrum — are necessary to explore more 
thoroughly the dual cascade of energy to small scales 
in both physical and velocity space. As discussed in 
Sec III HI this inherently kinetic and nonlinear mechanism 
plays a key role in the ultimate conversion of turbulent 
fluctuation energy into heat. 



V. APPLICATIONS 

Inertial range turbulence arises across broad range of 
space, astrophysical, and laboratory environments. The 



dissipation of the turbulence generally occurs at scales 
smaller than the ion mean free path; therefore, kinetic 
plasma physics determines the inevitable plasma heating 
due to the dissipation of the turbulence. In this section, I 
will identify a number of systems in which an understand- 
ing of the kinetic cascade is necessary to determine the 
dynamics responsible for the flow of energy and plasma 
heating. Identification of the characteristic length scales 
for a specific system — including the driving scale L, the 
ion mean free path A^^-^^p^, and ion Larmor radius pi — 
helps to locate its place within the theoretical framework 
of the kinetic cascade presented in Sec [TTl 

The most accessible turbulent astrophysical plasma is 
the solar wind; the ability to make in situ satellite mea- 
surements makes this an invaluable environment for the 
study of the kinetic cascade. Measurements of turbulence 
in the solar wind show a —5/3 magnetic energy spectrum, 
interpreted as the beginning of a turbulent inertial range, 
at scales L ~ 10" to lO^^ cm [H [H, HO, [Mi . The so- 
lar wind is very weakly collisional, with a mean free path 
of approximately 1 AU, or ^^aipi ^ ^^^^ '^'^ l63j . 
The ion Larmor radius in the solar wind is pi ~ 10^ to 
10^ cm [sl, i6J], yielding pi/L ^ 10^^ and giving an in- 
ertial range over five orders of magnitude in scale. This 
rather small inertial range, relative most astrophysical 
contexts, means the threshold at which cj — > il^ must 
be carefully evaluated; a more thorough and quantitative 
discussion of the applicability of gyrokinetic theory to the 
turbulent solar wind is contained in Howes et al. [53| . Al- 
though driven at a coUisionless scale L < Aj^j-p^, one ex- 
pects the general picture of the kinetic cascade described 
here to apply, as discussed in Sec III II It is important to 
note, however, that the coUisionality is sufficiently weak 
in the solar wind that the equilibrium particle distribu- 
tion functions are observed to deviate significantly from 
Maxwellian[6^ [63j . Therefore, velocity-space anisotropy 
instabilities, such as the mirror and the firehose, may 
play a role in the transfer of energy outside the scop e of 
this paper, as suggested by satellite observations [651 [66j. 

Interstellar scintillation due to electron density fluc- 
tuations in the interstellar medium (ISM) of our galaxy 
demonstrate a power law behavior consistent with an ex- 
ponent of —5/3 over twelve orders of magnitude in scale, 
from L 10^" to 10* cm These electron density 

fluctuations are considered to be passively mixed by the 
Alfven wave turbulence in the plasma [ij, [H, [s^l , and 
thus can be used as a probe of the plasma turbulence. 
The low end of the range is consistent with the ion Lar- 
mor radius in the ISM, pi ^ 10* cm The mean free 
path in the ISM is estimated to be ^Ynipi ^ ^^^^ 
The turbulent ISM therefore ideally fits the picture of 
inertial range turbulence in kinetic plasmas described in 
this work. 

In the accretion disk around a black hole, the mag- 
netorotational instability (MRI) [g^ is able to tap the 
energy of the gravitational potential to drive turbulence; 
this turbulence leads to enhanced angular momentum 
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transport, allowing the plasma to accrete onto the black 
hole. The MRI injects energy into the turbulent cascade 
on a scale comparable to the scale height of the thick ac- 
cretion disk, L ^ 10^'^ cm [70*1; this turbulence cascades 
down to the scale of the ion Larmor radius pi ^ 10^ cm. 
Due to the predicted high temperatures and low den- 
sities of the accretion disk, the mean free path in the 
plasma is of order the system size, ^^aipi ^ ^^^^ cm. 
The radiation emitted from the hot, magnetized plasma 
depends on the black hole properties and the heating of 
the plasma species by the kinetic dissipation of the tur- 
bulence [tiI, [zl, [zl] . To understand X-ray observational 
data from Chandra, the thermalization of the turbulent 
energy by the kinetic cascade must be characterized. 

Although the gradient-driven turbulence studied in 
the context of fusion plasmas does not develop an iner- 
tial range, energetic ions in a burning fusion plasma — 
arising either from Neutral Beam Injection or as hot 
a-particles generated by the fusion reaction — can drive 
shear Alfven wave instabilities at scales substantially 
larger than the ion Larmor radius of the thermal plasma 
74 [75I . ItgI . IttI . ItsI I . Experimental measurements have 
confirmed the generation of shear Alfven modes both by 
energetic ions |79l . [80| and by thermal ions [sTl . The 
wavelength of the Alfven modes driven by energetic parti- 
cles is comparable to the size of present-day experiments, 
so a discrete few Alfven eigenmodes arise; for the larger 
scale plasma of the International Tokamak Experimental 
Reactor (ITER), a broad range of these modes may be 
driven, possibly leading to the development of a turbu- 
lent cascade. Although the primary focus of early stud- 
ies has been on the confinement of the energetic ions in 
the presence of these Alfven wave instabilities, the ki- 
netic cascade driven by energetic ions may prove to be 
an important channel by which energy from the fusion a- 
particles is transferred to the kinetic scale of the thermal 
plasma and converted into heat. 

VI. CONCLUSIONS 



tion 

2. Collisions are required to increase the entropy and 
thus achieve irreversible heating of the plasma. 

3. The ultimate thermalization of the turbulent en- 
ergy is achieved in weakly coUisional plasmas 
through an entropy cascade, a dual cascade of en- 
ergy to small scales in both physical and velocity 
space. 

I have outlined the numerical approach chosen to in- 
vestigate turbulence in the transition regime from the 
Alfven wave to the kinetic Alfven wave cascade. The 
computationally demanding simulations use a Langevin 
antenna to model driving from larger scales in the in- 
ertial range and a hypercollisional operator to remove 
energy at the smallest resolved scales. Results from a 
simulation with high velocity-space resolution are con- 
sistent with the development of v±_ structure due to the 
action of the entropy cascade. 

Inertial range turbulence arises in kinetic plas- 
mas across a broad range of space and astrophysi- 
cal plasmas — including the solar wind, the interstellar 
medium, and accretion disks around compact objects — 
and may play an important role in the thermalization of 
fusion a-particle energy in next-generation burning plas- 
mas. 

A program of nonlinear gyrokinetic simulations of ki- 
netic, inertial range turbulence is necessary to explore the 
key elements of the kinetic cascade described here. The 
rich plasma physics involved in the inherently kinetic and 
nonlinear nature of the predicted entropy cascade, along 
with the computational intensity of the numerical work, 
make this important physics problem both fun and chal- 
lenging. Additional aspects of kinetic turbulence must 
also be examined, such as the effect of strong fast wave 
turbulence on the Alfven, slow, and entropy modes; this 
requires a treatment beyond the scope of gyrokinetics. 



In this paper, I have presented a general picture of a 
fundamental plasma physics process: the transfer of tur- 
bulent energy through an inertial range from the driving 
scale to dissipative scales followed by the conversion of 
this energy into heat. The details of the kinetic cascade of 
the generalized energy are not yet well understood. Nu- 
merical simulations of kinetic turbulence will play a cru- 
cial role in illuminating the key elements of the cascade. 
Here I emphasize several of these important elements: 

1. The nonlinear cascade of energy in kinetic plas- 
mas preserves a generalized energy given by equa- 
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